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Abstract 

This article is devoted to the study of two-dimensional Bose gases harmonically confined. We first summarize 
their equilibrium properties. For such a gas above the critical temperature, we also derive the frequencies and the 
damping of the collective oscillations and we investigate its expansion after releasing of the trap. The method is well 
suited to study the collisional effects taking place in the system and in particular to discuss the crossover between 
the hydrodynamic and the collisionless regimes. We establish the link between the relaxation times relevant for 
the damping of the collective oscillations and for the time-of-flight expansion. We also evaluate the collision rate 
and its relationship with the relaxation time. 
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Version frangaise abregee Cet article est consacre l'etude des gaz de Bose a deux dimensions en presence d'un 
confinement harmonique. Nous abordons tout d'abord les proprietes d'equilibre de ce systeme. Nous derivons 
ensuite, au-dessus de la temperature critique, l'expression des frequences et de l'amortissement des modes collectifs 
de basse energie et nous etudions avec le meme formalisme revolution du nuage d'atomes lorsque le confinement 
est brutalement supprime. La methode utilisee permet de decrire le gaz dans tous les regimes collisionels, du regime 
sans collision au regime hydrodynamique. Nous etablissons le lien entre les temps de relaxation qui decrivent les 
modes d'oscillation et l'expansion du nuage apres coupure du piege. Nous evaluons egalement l'expression du taux 
de collision et sa relation avec le taux de relaxation des modes de basse energie. 
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1. Introduction 



Quantum gases in reduced dimensionality are now experimentally available. Two-dimensional Bose gases have 
been realized by trapping atomic hydrogen at the surface of liquid helium [1,2]. Such a gas, confined into a box, 
undergoes at sufficiently low temperature a superfluid transition known to be of the Berezinskii-Kosterlitz-Thouless 
type [3,4,5]. 

Recently, a few experiments with laser-cooled atoms have approached the two-dimensional regime. The method 
used consists in realizing very anisotropic confinement in such a way that one degree of freedom is frozen to 
zero motion oscillations. The crossover to two-dimensions occurs when the thermal energy is below the vibra- 
tional energy in the tightly confined direction. Such a regime has been reached in standing-wave dipole traps 
[6]. Alternatively, an elliptically focused laser beam has been fed by a Bose- Einstein condensate leading to a 
two-dimensional configuration [7]. Evaporative cooling of an atomic gas has also been performed at the crossover 
to two dimensions in an optical surface trap [8]. A two-dimensional Bose- Einstein condensate has been achieved 
with ultracold cesium atoms trapped in a gravito-optical surface trap [9]. Notice that two-dimensional trapping 
by means of a field-induced adiabatic potential has been recently proposed [10] and realized experimentally in the 
group of H. Perrin and V. Lorent [11]. In all those systems the two-dimensional confinement can be considered as 
harmonic. 

The harmonic confinement introduces very different features from the physics one would obtain in a box. As 
explained below, for an ideal Bose gas, Bose-Einstein condensation occurs at finite temperature in the presence 
of a harmonic potential in contrast to the homogeneous case [12]. 

The paper is devoted to the physics of two-dimensional Bose gases harmonically trapped above the critical 
temperature. First, we discuss the thermodynamic properties and the shift of the critical temperature of an ideal 
gas at fixed number of particles due to the finite size of the sample. We discuss briefly the role of interactions. 
Second, we study the dynamics of the gas through low- lying collective modes and time of flight expansion by 
including dissipative and mean-field effects. 



2. Thermodynamical properties 
2.1. Ideal gas 

We consider N bosons confined by a two-dimensional harmonic trap of angular frequencies lo x and u> y . By 
contrast with a non-confined Bose gas, Bose-Einstein condensation is expected to occur since the number of 
atoms in the excited states saturates to an upper value N^^: 

N > 1 ^ 7T 2 k B T k B T 

max ^—f exp[h(n x u) x +n y LUy)/(k B T)] — 1 ~~ 6 hcj x hu y ' 

The critical temperature for an ideal Bose-Einstein condensation T° is readily obtained from this formula by 
setting iV^ax = N: 

with (D = {iJxLUy) 1 / 2 . Strictly speaking, this value for the critical temperature is valid only in the thermodynamic 
limit (a) — > 0, N — > oo with the product u) 2 N keeping a constant value). We stress that due to the confining 
potential, "true condensation" occurs in the sense that the first excited states have a probability of occupation 
which tends to zero in the thermodynamical limit. 

We find worthwhile to examine the finite-size correction [13] to the critical temperature as real experiments 
are more likely carried out with a small number of particles. The scaling of the expected correction has been 
investigated in [14], we provide in the following a more quantitative estimate. For a finite number N of particles, 
the "critical temperature" T c , still defined as the temperature for which iV^ax = N, is shifted downwards. One 
can calculate analytically, for an isotropic harmonic confinement, the leading term correction for large but finite 
N by evaluating the sum (1) more accurately (see appendix A): 
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Figure 1. Critical temperature T c (normalized to T®) for a finite number N of atoms. The solid line is the analytical result given 
by (3), while the circles arc given by numerically solving for T the equation N max (T) — N. The agreement is excellent even for low 
values of N. 



The prediction of Eq. (3) is compared, on Fig. 2.1, to the exact value of T c obtained by solving numerically 
the equation iVmax(r) = N. We notice that the shift in T c is quite small by contrast with its three-dimensional 
counterpart, even for N as low as a few hundred. 

2.2. Interactions 

The interatomic collisions in a tightly confined Bose gas have been studied in detail in [15]. We summarize in 
this section their main results. Similarly to the three-dimensional scattering problem, the scattering properties 
in two dimensions can be formulated by means of a single length d* . This characteristic length depends on the 
detailed shape of the interatomic potential. For realistic momenta q of particles in ultracold gases we always have 
the inequality qd* <C 1. The expression for the two-dimensional differential cross section for identical bosons reads: 
da/d6(q) = \h(q)\ 2 /(Anq), where q = m(\r 2 — Vi)/(2h) is the relative wavevector of the colliding particles and the 
dimensionless function h(q) depends logarithmically on q. The mean field potential in the non degenerate regime 
is U(r) = 2gn(v) with n the density, g = h 2 h{q)/m and q = (mfc_BT) 1//2 /fi the mean thermal wave vector. 

The departure from the ideal gas behavior is measured through the ratio between the mean field energy and the 
level spacing. For weak interactions, namely gn <C to, interactions can be taken into account perturbatively and 
the physics of the ideal gas is valid. In the opposite limit gn to, interactions brings about major differences. 
One enters a regime of quasi-condensate [16]. From thermodynamics, one readily establishes the implicit equation 
for the density: 



where A = h(2ivmkBT)~ 1 '' 2 is the thermal de Broglie wavelength, /3 = 1/fcsT and L/trap is the trapping potential 
energy. Actually, in the non-perturbative interactions regime, the local density approximation is valid. Repulsive 
interactions tend to push the gas towards the external region, thereby reducing the density of the gas. One expects 
that this effect causes a natural decrease of the critical temperature. This is indeed the case and in contrast with 
its three dimensional counterpart, Eq. (4) is soluble at all temperatures [17]. To discriminate which state with 
or without condensate the system will prefer, the authors of Ref. [18] compare the free energy in both situations 
and predict that low-energy phonons destabilize the two-dimensional condensate. The nature of the phase in 




(3) 



„( r ) = --Llog(l - e-^WrHM-^ 



(4) 
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the degenerate regime remains somewhat an open question [19,16,18]. One may wonder if this state involves 
superfluidity at sufficently low temperature. 



2.3. Expansion of the density 



In on-going experiments on Bose-Einstein condensates of alkali atoms, the density profile is readily obtained 
by in situ images of the cloud. We give in this section the corrections to the density up to the second order due 
to the mean field contribution. However, this expansion is valid not too close to degeneracy. We expect that the 
repulsive mean field contributes to the decrease of the density. We capture this effect perturbatively by performing 
an expansion of the density profile with £ = 2gf3\~ 2 as the small parameter. £ is the ratio between the mean field 
energy close to degeneracy and the temperature. By denoting the equilibrium density in absence of mean 
field, we work out the expansion: 

n (r) = 4 0) (r) +Cn ( ( , 1) (r) +C 2 n< 2) (r) + ... 
n<j 0) (r) =-A- 2 log(l- K ) 

n W (r)= _n^_^_ l lQg(1 _ K) y (5) 
with k = «(r) = exp(— /3(Z7 tr a P (r) — /x)). 



3. Dynamics 

3.1. Quantum Boltzmann equation 

The dynamics of the gas is described by a two-dimensional quantum Boltzmann equation: 

§£+v 1 -V r / + ^-V vl / = / coI1 [/] (6) 
where / is the phase space distribution function. The explicit form for the collision term is 

I coll [/(r, Vl )] ~ ^ l -^ff- J d6 d 2 v 2 [/(r, Vl ,)/(r, v 2 ,)(l + /(r, vi))(l + /(r, v 2 )) 
-/(r, vi)/(r, va )(l + /(r, Vl /))(1 + /(r, v a /))] 

It accounts for elastic collisions between bosonic particles 1 and 2, with initial velocities vi and V2, and final 
velocities Vj./ and v 2 <, being the scattering angle in the center of mass frame. In the next sections, we first derive 
the expression for the collision rate with the interplay between statistics and mean field. Then, we investigate 
the role of collisions (mean field and relaxation) on the frequencies of the low-lying modes for a two-dimensional 
harmonic trap and we derive the equations for expansion after an abrupt switching off of the trap. The method 
used to solve those latter problems relies on a scaling ansatz. It is first introduced in the context of the collision- 
dominated hydrodynamic regime, and then adapted to the quantum Boltzmann equation [21]. 



3.2. Collision rate 



The collision rate plays a crucial role for the relaxation dynamics, it permits to compare the mean free path to the 
size of the cloud and to characterize the collisional regime of the sample between collisionless and hydrodynamic. 
Well above the critical temperature, the mean field and the bosonic statistics can be neglected and the equilibrium 
distribution reads [20]: 

/o(r ,v)=iV.^.^ ; e-^(r.v)ABT (7) 
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with 



H (r,v) = ^rnv 2 + ^mw^ 2 + imcj 2 i/ 2 . 



(8) 



We readily obtain the expression for the collision rate 70 in this limit by integrating the first term of the collision 
integral of the classical Boltzmann equation for the equilibrium function fo and dividing by the number of particles: 



7o = 



m 3 |fe(g) | s 
h 3 2-kN 



J d 2 rd 2 Vl d 2 v 2 Mr.vO/ofrv,,) = ^^\h(q)\ 2 - 



(9) 



Closer to degeneracy but still above the critical temperature, one expects that the result is modified by two 
opposite effects. On the one hand, the mean field tends to decrease the density. On the other hand, the statistics 
tends to shrink the cloud. To evaluate the competition between those effects, we restrict ourselves in the following 
to the case of an isotropic harmonic confinement for sake of simplicity. We first expand the phase space distribution 
function: /(r, v) = /o(r, v)[l + / (r, v) — 2(3gn (r)]. From the quantum Boltzmann equation, the collision rate 
can be expanded to the first order 



7 = 



m a \h(g)\ 2 
h s 4n 2 N 



J d6 d 2 v 2 d 2 v 1 d 2 r f (r, vi)/o(r, v 2 ) 



1 + 



/o(r, vi) + / (r, v 2 ) + /o(r, v v ) + / (r, v^) - 4g/3n (r) 



(10) 



The signs in (10) clearly show the enhancement of the collision rate due to the statistics and its decrease due to the 
repulsive mean field. After a lengthy but straightforward calculation, we have worked out the explicit expression 
for the expansion of the collision rate with respect to statistics and mean field contribution: 
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7 = 7o 



1 + ^1- 



gm_ 

irh 2 



Actually, the collision rate cannot be measured directly. However, all dissipative dynamics involve the collisions. 
We propose in the following to extract information on the collisional regime through the study of the collective 
oscillations or of the time-of-fiight expansion, which can be both investigated experimentally. 



3.3. Scaling ansatz for the hydrodynamic regime 



In this section, we determine the frequencies of the collective oscillations of an harmonically trapped two- 
dimensional gas, and the equations for its time of flight when the confinement is switched off. For this purpose, 
we transpose the hydrodynamic approach developped in [22] for a three dimensional trapped Bose gas to a two- 
dimensional one. The scaling ansatz for the density reads: n(r,t) = no(R(t))/f2o(£), where no is the equilibrium 
density, R = A(t).r with r the coordinate and A(t) a diagonal matrix with time dependent scaling coefficient 
A — di&g(l/a x (t),l/a y (t)). The parameter cn gives the dilation along the ith direction. One readily obtains 
flo{t) = (det(A))~ 1 from the normalization J d 2 rn(r,t) = J d 2 mo(r). Note that the scaling solution makes sense 
only if A has non vanishing coefficients. The local hydrodynamic velocity field v loc is obtained from the ansatz for 
the density through the equation of continuity, which yields: v loc = — A~ 1 A.r, where A — dA/dt. The equations 
for the coefficients of the matrix A are obtained from the Euler equation. For collective oscillations, we find the 
following set of nonlinear equations: 

2 2 

2 LOx 2 

a x + uj x a x — = and a y + uj y a y ^— = 0. 

^x^y Oj y Q/x 

For a small deviation from equilibrium, the linearization of the system leads to two hydrodynamic frequencies of 
oscillation: 

3 + 3/3 2 =b (9 - 14/3 2 + 9/3 4 ) 1/2 

with p = uiy/ujx- For an isotropic trap (oj x = uj v = u>o), we find V^loo for the quadrupole mode and 2a>o for 
the monopole mode. For non isotropic trap, the two modes correspond to a superposition of the quadrupole and 
monopole modes. The equations for the time of flight in two dimensions are readily obtained in the same way: 
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a x = and a y = — — . 

This set of equation shows that a high initial collision rate with respect to the trap frequencies (implicitly assumed 
for the validity of hydrodynamic equations) leads to an asymptotic anisotropic expansion [21]. 

3.4. Scaling ansatz for the Quantum Boltzmann Equation 

In most cases the hydrodynamic formalism is not satisfactory for describing a Bose gas. Indeed, a time of flight 
expansion is accompanied by a dilution of the sample. As a result, the collision rate decreases. Hydrodynamic 
equations are mostly valid at the beginning of an expansion but certainly not after a long time [23]. We emphasize 
that the hydrodynamic regime is quite difficult to reach experimentally [24,25] since a high collision rate means a 
high density for which the inelastic collision rate is magnified. However by exploiting Feshbach resonance, it has 
been possible recently to enter the hydrodynamic regime for fermionic gases in a three dimensional trap [26]. In 
order to take into account the evolution of the density one needs to solve at least approximately the Boltzmann 
equation. This equation permits us to describe the crossover between the collisionless and the hydrodynamic 
regime. 

Following [21], we make the following ansatz for the distribution function: /(r, v,£) = /o(R, V)/fi(t) with 
fo the equilibrium distribution, V = B.(v — v loc ), B being a diagonal matrix B = diag(l/&i (i),l/b y (t)), 
and Q(t) = det~ 1 (AB). The parameter bi gives the effective temperature in the ith direction and fo obeys the 
equilibrium Boltzmann equation: 

V- f v dfo 2 „ df 2g dn df \ 

where n = (m 2 /h 2 ) J d 2 Vfo is the equilibrium density and go the strength of the interaction for the equilibrium 
temperature To. The scaling ansatz method does not provide an exact solution of the Boltzmann equation. 
However, it is in reasonable agreement with the numerical simulations based on molecular dynamics [21,27] and 
with experiments [28] . By substituting the scaling ansatz for the distribution function in the Boltzmann equation 
and taking into account the properties of the equilibrium distribution Eq. (11) one obtains: 



-y 



9o cnb 1/2 Hi a o I 9Ri \ 2b % ai I dVi 



2 g 1 \ Ri dfo 

ai+LOiai - 1 



go a z n, o-j J b 1/2 dVi 



/con + tJ/o- (12) 



This way, we restrict further the solutions. Such an approach would fail for instance to take into account rota- 
tion in presence of the mean field. Eq. (12) gives constraints on the scaling coefficients a* and h. By denoting 
(C)o = (m 2 /h 2 ) f C/o(R, V)d 2 Rd 2 V/N, we readily derive the following set of equations by multiplying Eq. (12) 
respectively by C — RiVi and C = V 2 and by performing the integration over the phase space: 

bl 9 1 1 )(V 2 )o + (&i + u?a,i - —-TT- - I - h ' ' - 



<'• II"/ \ '' II.". 

- j d 2 Rd 2 VV 2 I coll - (13) 



bi hi _ m 2 ft 

b~ + ~ l ^~h 2 ~Wy 2 )o 



In order to capture the physics of the collision integral, we treat this term within the relaxation approximation 
[21]. It allows us to recast the l.h.s of the last equation of (13) in the form: 

^ f j2 D ,2 1/1/ 2 r „ _ 1 I bi- bj 



dTRdrVVChow = -— ' J with i + j, 



h 2 N{V 2 )o J bi \ 2r 

where r is the relaxation time which corresponds to the average time between collisions. Using the properties 
of the equilibrium distribution, we finally obtain the following closed set of non linear equations for the scaling 
parameters: 
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... 2 2 h 2 1 / h g 1 
ai + u>; a, - hw t q — 



a, I a t g a % \[ i a,j 

b i + 2^6 i = -^(6i-6 J ) with j, (14) 

a» zr 

with £ = go(no)o/('mujf(Ri}o). One recovers the collisionless regime by taking r = co. In this limit, we have 
a simple relation between fei and at: hi = a~ 2 . In the opposite limit (hydrodynamic regime), local equilibrium 
is always ensured because ol the high collision rate. As a consequence, the contribution ol the collision integral 
vanishes because of local equilibrium and b x = b y = l/(a x a y ). As the mean temperature is constant g = go, we 
recover in this limit the set of equations (11). In the next sections, we use this set of equations to derive the 
frequencies of the collective oscillations and we establish the equations for a time of flight expansion. 



3.5. Collective oscillations 



Low lying collective oscillations are reproduced by the time dependent dilatation parameters a x and a y of the 
trapped cloud around equilibrium. As a consequence, the temperature dependence of the interaction strength can 
be neglected g = go. To evaluate the relaxation time to, at least at equilibrium and in the high temperature limit, 
we perform a gaussian ansatz in the collision integral as explained in details in [27] for a three dimensional system. 
We find to = 2/"fo- Expanding Eqs. (14) around equilibrium (o» = 6» = 1) we get a linear closed set of equations 
which can be solved by searching for solutions of the type e 1 "*. The associated determinant yields the dispersion 
law: 



A[uj] - —B[lo] = 0, 

T() 



(15) 



where A[uj] = u (lo — u c1+ )(lo ~uj c1 _), B[uj] = u(lo — LO hd+ )(io —u hd _) and (cl) and (hd) refer to the collisionless 
and hydrodynamic regimes respectively. The values for u>^± and cu^ d ± are given by 



2 Wi 



^hd± - T 



(1 + /3 2 )(4 - ± {(1 + /3 2 ) 2 (4 - 2 - 32/3 2 (2 - £)} 1/2 
3(1 + /3 2 ) ± {9(1 + f3 4 ) - 14/3 2 } 1/2 



(16) 



with (5 = u) y /u) x . 

The roots of A that correspond to the collisionless regime in presence of mean field coincide with the one derived 
by [29]. The frequency of the monopole mode [29] for an isotropic trap (lo x = uj y = luq) is found to be 2u>o whatever 
the collision regime. In two dimensions the collisions (mean field and relaxation) do not affect the frequency of 
this mode, this comes out from the fact that it is in the kernel of the integral of collision [27] . The symmetries 
of a two dimensional system for a contact interaction can explain this surprising result [30]. The hydrodynamic 
frequencies do not depend explicitly on the mean field, it is also a specific result of two dimensions low lying 
oscillations. As a consequence, the effect of the mean field upon collective oscillations is maximum in the crossover 
between collisionless and hydrodynamic regime. 



3.6. Time of flight 

In the following we investigate the evolution of the cloud in a time of flight. In this technique, the asymmetric 
trapping potential is switched off and the evolution of the spatial density is monitored. After a long time expansion, 
the aspect ratio gives information on the characteristics of the gas. Below the critical temperature, the anisotropic 
expansion is the favored signature of Bose- Einstein condensation. Actually, above the critical temperature and in 
the collisionless regime one expects an isotropic asymptotic expansion of the cloud reflecting the isotropy of the 
initial velocity distribution. However, anisotropic expansion also arises above the critical temperature, when the 
mean free path is small compared to the size of the cloud. 

An analytic approach has been proposed in the full hydrodynamic regime [22] for a three dimensional Bose gas 
and has been adapted to a two dimensional gas in section 3.3. Alternatively, the expansion of an interacting Bose 
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gas initially confined in a three dimensional trap above the critical temperature has been investigated by means 
of Monte Carlo simulations [31]. Two recent papers [21,23] propose a way to provide an analytical interpolation 
between the two opposite collisionless and hydrodynamic regimes. The authors of Ref. [23] divide the expansion 
in two stages: the first one is considered as fully hydrodynamic and the second one as collisionless. The authors 
of Ref. [21] use a scaling ansatz to solve approximately the Boltzmann equation. We adapt in the following this 
latter method to a two dimensional gas. 

The same procedure as the one used previously for the collective oscillations leads to the following set of 
equations: 

2bi , 2 J bi g(b) 1 _l\ 
k +2^-h = -—^—{b i -b j ) with i^j, (17) 

d; 2,T(CLi) 

where b = (b x + b y )/2. The confinement does not appear explicitly as in (14). The dependence with the parameters 
di and bi of the collision term through the interaction strength g and the time relaxation r has to be taken 
into account. As recalled in section 2.2, g depends logarithmically on the mean temperature T = Tob. Since 
the collision rate scales as the density, we deduce the scaling dependence of the relaxation time: r = To\\ i a i . 
During the expansion the scaling parameters bi decrease and as a consequence the gas cools down, revealing the 
isentropy of the expansion. By contrast with three-dimensional expansion, the mean field contribution is slightly 
magnified during the expansion through the dependence of g upon the mean temperature. The persistence of the 
hydrodynamic expansion that may happen initially is more pronounced with respect to the three-dimensional case 
because of the more favorable dependence of the relaxation time on the scaling parameters. 



4. Conclusion 



In conclusion, our work shows that the coherent part of collisions included through a mean field term can 
manifest itself in various experimentally measurable ways already above the critical temperature: by modifying 
the frequencies of low-lying modes, or in competition with dissipative effects by producing an anisotropic expansion 
for an initially anisotropic trapping. 



ACKNOWLEDGMENTS We acknowledge fruitful discussions with M. Holzmann and useful comments from 
H. Perrin. 



Appendix A. Finite size effect on the critical temperature T c 

For N bosons confined in a two-dimensional isotropic harmonic oscillator of angular frequency to, the critical 
temperature T c is the solution of the following equation 

N _ ^ n + 1 

~ ^ exp[fio;n/(fc B T c )] - 1 ' ( ' ' 

Since T c is close to T c ° = SwVW/(ffe), one has e = hu/(k B T c ) ~ N~ 1/2 < 1. 

We therefore expand Eq. (A.l) for £<1, retaining only the lowest order terms, up to 1/e. The sum can be 
replaced by an integral, provided we add a remainder TZ(e) to keep an exact expression: 

oo 

N -J^h dx+n{£) - (A - 2) 
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The expansion of the integral is straightforward 

oo 

X + 1 -d*=Zfi-i-^ + 0(l) (A.3) 



exp(ex) — 1 6 e 2 e e 

1 

and the remainder can be evaluated with the Euler-McLaurin asymptotic formula [32] 

7e( £ ) = -gt/^(D with /(x) = ^J_ (A . 4) 

where -Bj, is the fcth Bernoulli number. One then expands the derivatives of / up to order 1/e, and one finds: 

ft( e ) = | +0(1), (A.5) 

where the coefficient S is given by the following series: 

S = 1 + 1/12 -1/120 + 1/252 + ••• (A.6) 

which turns out to be a diverging series. Nevertheless, retaining only a finite number of terms of this series actually 
gives a very good approximation [33] of the value of S ~ 1.079 ± 0.003. Therefore, 

., 7T 2 1 0.079 lne 

7V= T - + ^- — + 0(1). (A.7) 



Substituting e by ty 2 /6N (1 — 5T C /T^ ) in the previous equation, and expanding up to first order in <5T C /T° finally 
yields the shift of critical temperature: 

ST C 0.195 In N- 0.066 

n~ 7W ' (A - 8) 
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